Phantom fields: bounce solutions in the early Universe and S-branes 
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The cosmological model with two phantom scalar fields with the special choice of field's potential 
is considered. The obtained regular solution describes a bounce with a subsequent transition to the 
de Sitter stage of the expansion of the Universe. This solution also could be interpreted as a S-brane 
. solution. 
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C — ■ ' Recent observations of Type la supernovae show an accelerating Universe at the present epoch [l|, 0] . Such acceler- 
ated expansion is caused by presence of some antigravitating substance which act as repulsive force in the Universe. 
The immediate observation of this substance is sufficiently hard due to its weak interaction with electromagnetic 
radiation. In this connection this substance was called dark energy (DE). The true nature of DE remains a puzzle 
till now. The simplest model is a description of DE as a constant A-term (the so-called ACDM model, see [2j, for 
^v^j ' example). However, at the present time it is rather difficult to find physically proved explanation of a value and 
(T^) , nature of the A-term. 

In this connection, consideration of dynamical models of DE with the varying equation of state seems to be more 
natural. In this case the effective equation of state of DE corresponds to w < —1/3, while for the ACDM model 
w = — 1. All the models laying in these limits on w can be divided into two big categories: 1) the source of DE 
is different scalar fields (quintessence models, Chaplygin gas, etc. [3, 0, 0]); 2) presence of DE is a consequence of 
O modification of the Einstein gravitational equations by addition of terms depending on scalar curvature, invariants of 
the Ricci and Riemann tensors, etc. @, 0, [1, @- Consideration of the evolution of the present Universe within the 
framework of these models allows to achieve the certain success in quantitative comparison of theoretical results with 
observations. 

However, the recent observations [Iol | have shown that w < — 1 (see Ref. fill ] where this phenomenon is described 
by a model- independent way). Such a type of DE, also known as phantom energy can be introduced within the 
framework of both directions mentioned above. In particular, in case of scalar fields phantom energy is realized by 
change of a sign before a kinetic term in a field Lagrangian. Choosing, further, a certain form of potential energy of 
this field, it is possible to get a satisfactory agreement with observations by selecting parameters of a model (see, e.g., 
[l3| and references therein). On the other hand, one can try to model an inflation in the early Universe using phantom 
fields [IH and consider quantum evolution of the early Universe within the framework of phantom cosmology [loT ]. 

Phantom DE for two scalars (where only one of scalars is phantom) was considered in [16| . The idea of unification of 
inflation with acceleration for phantoms was proposed in [l7| • The paper [l8[ is devoted to the analysis of conformally 
coupled phantom scalar fields. The authors apply the Ziglin theory of integrability and find that the flat model is 
nonintegrable. It means that we cannot expect to determine simple analytical solutions of the Einstein equations. 
The metric presented here shows the same: a solution of the Einstein equations coupled with two interacting scalar 
fields can be found only in a numerical form. General review on modified gravity as DE one can find in 

Another possible application of phantom scalar fields consists in modelling of the so-called spacelike branes (S- 
branes) with their help. The S-branes are topological defects localized on a spacelike hypersurface. An example of 
such a solution is a domain wall in a 4D scalar field theory [2(|. In this model the theory of one scalar field with 
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effective potentiai energy with spontaneous symmetry breaking is under consideration. This effective potential has 
two minima. For a static case, this problem has the well-known solutions in the form of kinks [2lJ . These solutions 
start from one minimum (at x — > — oo) and tend asymptotically to another minimum of the effective potential (at 
x — > +00). However, if one consider phantom helds, the similar solutions in a form of kink will be obtained not for 
the space variable x but for the time t. This so-called S2-brane solution is a timelike kink [22| . 

In this paper we consider a model of S-brane formed by two phantom fields having local and global minima. Below 
we show that solutions similar to the soliton-like solutions mentioned above could exist and for a case when an effective 
potential is choosing in the form of two interacting scalar fields ([3]) . This potential is an analog of the well-known 
model with one scalar field with a Mexican hat potential. The fundamental difference of these two potentials consists 
in presence of two scalar fields in ^ that provide existence not only kink-like solutions but also another types of 
soliton-like solutions. 



II. EQUATIONS 

Let us consider the flat cosmological Friedmann model with the metric: 

ds 2 = dt 2 - a(t) 2 (dx 2 + dy 2 + dz 2 ) 
The Lagrangian for phantom scalar fields <p, x plus gravitation is: 

R 



where 



L = 



L( ._ r , - ^w-^x-n^x), 



(1) 



(2) 



(3) 



Here Vq is defining by initial conditions and can be considered as a cosmological constant A. Introducing the dimcn- 
sionless curvature tensor p\ — l 2 iR\, the Einstein and scalar field equations can be written as 
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(4) 



(5) 
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dx 



(6) 



where G is the gravitational constant, l p i - the Planck length. Using the metric {T]) and introducing dimensionless 
variables and constants <f> — > (fif^/Swlpi, x ~^ x/v / 87r^/, t — > \f^nl v it^ mi — > mi/^/8Trl p i, 777,2 - ► Tn^/V^lpi, A1.2 — > 
2Ai,2 one can rewrite the equations (J3j)-(|6j) as (in the Planck units, i.e. at c = H = G = 1): 
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~tf-x 2 + v{^x) 



(7) 



(8) 



■3^ = ^[^ + Al (0 2 - 777?)], 

a 



X + 3-x - X W + A 2 (x 2 " m 2 )} , 



(9) 
(10) 
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where 



and 



7(0, X) = y (0 2 - m\f + yV ~ m\f + 0V - V , (11) 



v = y^o- m i) 2 + t (xS " mlf + ^ (12) 

The last value follows from Eqs. ([7]), (fTT|) and the initial conditions (fT3ll - (fT5|) . (The equation ((8]) was obtained by 
subtraction of the (I) component of (HJ from the {%) component.) For these equations, there are following initial 
conditions: 

(13) 
(14) 
(15) 

In Ref. [HI the research of mathematically similar equations have been carried out. Those equations describe the 
thick brane in the Randall-Sundrum scenario. It was shown that at choice of the potential in the form (|11[) . there are 
the regular solutions describing the thick brane. Such solutions can be obtained only if the potential has local and 
global minima. In the context of Ref. [23| , the specified scalar fields could be interpreted as a condensate of a gauge 
field. 

Having a purpose to obtain regular cosmological solutions, we choose phantom fields with the given potential as a 
model of scalar field. In this case Eqs. (f7))- (|TU]) arc mathematically similar to the equations describing the thick brane 
in Ref. [H. 



a(0) = 


ao, 


6(0) = 


0, 


0(0) = 


00, 


0(0) = 


0, 


x(o) = 


Xo, 


x(o) = 


0. 



III. NUMERICAL ANALYSIS 



For the numerical calculations we choose the following values of the parameters and initial conditions: 

a = 1, 0o = V3, xo = VlA, Ai = 0.1, A 2 = 1.0. (16) 

We apply the methods of step by step approximation for finding of numerical solutions (the details of similar calcu- 
lations can be found in Ref. [23|). 

Step 1. On the first step we solve Eq. ([9]) (having zero approximations ao(t) — aoe*, Xo(*) = V0~l5/ cosh 2 (i/4)). 
The regular solution exists for a special value toJ i only. For m\ < m\ i the function <fii(t) — » +oo and for m\ > m\ i 
the function 0j(i) — > — oo (here the index i is the approximation number). One can say that in this case we solve a 
non-linear eigenvalue problem: 0*(i) is the eigenstate and to* i is the eigenvalue on this Step. 

Step 2. On the second step we solve Eq. (fT0|) using zero approximation ao(t) for the function a(t) and the first 
approximation (f>*(t) for the function 0(t) from the Step 1. For to 2 < m\ 1 the function xi(^) ~~ * +°° an d for to,2 > m\ 1 
the function xi(t) ~ * — oo. Again we have a non-linear eigenvalue problem for the function XiM an d rn\ 1 . 

Step 3. On the third step we repeat the first two steps that to have the good convergent sequence (p* (t) , \* (t) ■ 
Practically we have made three approximations. 

Step 4. On the next step we solve Eq. (j5J| which gives us the function a\(t). 

Step 5. On this step we repeat Steps 1-4 necessary number of times that to have the necessary accuracy of 
definition of the functions a*(t), (f>*(t), ~X*{t) and the eigenvalues ml 2 . 

After Step 5 we have the solution presented in Figs. H]), ©. These numerical calculations give us the eigenvalues 
to* = 2.98530175, to*. = 2.538974725 and eigenstates a*(t), 0* (i),X*(*) • 



IV. CONCLUSIONS 



The toy model of inflation in the early Universe with account of the phantom fields is considered. The possibility 
of obtaining of bounce solutions is shown in principle. In this case the model evolves from a stage of contraction 
(a < 0), passes through a minimum of the scale factor (ao ~ 1 in the Planck units) and appears on the de Sitter stage 
of expansion (Fig. ([2])). Such evolution is provided in the model by a special form of the potential for the phantom 
fields Q3J). 
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FIG. 1: The graphs for the scalar fields (f>*(t),x*(t)- Both 
functions tend asymptotically to the local minimum at <j>* (t) = 
mi,X*(t) = 0. 
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FIG. 2: The Hubble parameter H = a/a tends asymptotically 
to the de Sitter stages: if H < - the stage of contraction, if 
H > - the stage of expansion. 



From the mathematical point of view, the problem was reduced to finding of eigenfunctions of the nonlinear 
differential equations (f8|)- (|T0|) . Only at the obtained values of to*,TO2, the solutions are regular and a bounce and an 
inflation take place. At other values of these parameters, the solutions are singular. 

Let us note that at the obtained values of the parameters mi 2 = m\ 2 , the de Sitter stage occurs over an infinite 
period of time. But if \m\^ — m\ 2 | <C 1 than there is a finite time of inflation of the Universe. 

The obtained solutions could be interpreted as S2-brane solutions for a case of two interacting phantom scalar 
fields. But in contrast to the usual kink-like S2-brane solutions [20, 22], our solutions are another type of soliton-like 
solutions. If the kink solution starts from one minimum of a potential and tends asymptotically to another one, our 
solutions start and finish in the same minimum (in the case under consideration, in the local minimum </> = mi, x — 0, 
see Fig. ([]}). By the terminology of [2l|, the kink-like solutions refer to topological solutions and our solutions - to 
non-topological ones. In contrast to the kink-like solutions which exist only in a case of presence of two minima, the 
existence of non-topological solutions is also possible at one minimum only. 
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